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We analyze the properties of a Luttinger liquid under the influence of a periodic driving of the 
interaction strength. Irrespective of the details the driven system develops an instability due to a 
parametric resonance. For slow and fast driving, however, we identify intermediate long-lived meta- 
stable states at constant time-averaged internal energies. Due to the instability perturbations in 
the fermionic density are amplified exponentially leading to the buildup of a superlattice. The mo- 
mentum distribution develops a terrace structure due to scattering processes that can be associated 
with the absorption of quanta of the driving frequency. 
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I. INTRODUCTION 

Due to remarkable progress in experiments it is within 
the scope of present technology to implement and simu- 
late the dynamics of quantum many-body systems with 
a high degree of controllability on the system parameters 
even under nonequilibrium conditions. While Quantum 
Dots provide the framework for the study of quantum 
impurity systems 1 , ultra-cold atoms in optical lattices 
constitute the basic structure for the experimental real- 
ization of interacting quantum many-body systems such 
as the Bose- or Fermi-Hubbard model 2 . 

In principle, one can imagine a variety of differ- 
ent nonequilibrium driving protocols. By now the 
nonequilibrium dynamics following interaction quenches 
or ramps is studied in great detail 3 . The influence of a 
periodic driving, however, on strongly correlated many- 
particle systems poses new challenges especially in terms 
of methodology and has only been studied for a restricted 
class of systems. There has been considerable work on 
periodically driven strongly correlated impurity systems 
such as Anderson impurity 4 and Kondo models 5 ' 6 . Be- 
sides impurity systems, Falicov-Kimball 7 ' 8 and Hubbard 
models 9 ' 10 have been analyzed on the basis of a nonequi- 
librium extension of dynamical mean-field theory. Pe- 
riodically driven systems of interacting fermions in one 
dimensions have been studied for small system sizes 11 
and in the Luttinger liquid limit 12-14 . Periodically driven 
one-dimensional Bose-Hubbard models have been inves- 
tigated by exact diagonalization for small systems 15 and 
for large systems based on the time-dependent density 
matrix renormalization group 16 . 

In this work we investigate the nonequilibrium dynam- 
ics of interacting fermions in one dimension within a 
Luttinger liquid description induced by a periodic time- 
dependence of the interaction strength. The impact of 
a periodic modulation of the Fermi velocity onto the 
fermionic momentum distribution has been investigated 
recently 12 . Due to the above mentioned complexity of 
periodically driven systems it is instructive to explore 
those particular cases where exact and nonperturbative 



solutions are accessible such as in the case studied in this 
work. 

The periodically driven Luttinger liquid shows an in- 
stability in the long-time limit due to a parametric reso- 
nance 12 ~ 14 . However, on intermediate time scales meta- 
stable steady states of constant time-averaged energy 
densities can form with long lifetimes in case of slow and 
fast driving. In the adiabatic limit the instability occurs 
after a fixed number of periods independent of the driving 
frequency. Thus it is impossible to adiabatically follow 
the ground state of the system for a large number of peri- 
ods irrespective of how slowly the system is driven. The 
periodic driving promotes fermionic scattering processes 
under the absorption of quanta of the driving frequency 
yielding a terrace structure in the fermionic momentum 
distribution. In the slow driving limit the system be- 
comes unstable against perturbations in the fermionic 
particle density resulting in an exponential amplification 
of the perturbation in consequence of the parametric res- 
onance. 

This paper is organized as follows. In Sec. II we intro- 
duce the model of a periodically driven Luttinger liquid 
whose dynamics including the parametric instability we 
analyze in Sec. III. In Sec. IV we identify a meta-stable 
state for fast and slow driving based on the study of the 
internal energy density. The impact of the instability on 
the fermionic density and the momentum distribution is 
analyzed in Sec. V and Sec. VI, respectively. 



II. PERIODICALLY DRIVEN LUTTINGER 
LIQUID 

Consider a system of ID fermions of length L whose in- 
teraction strength is varied periodically with a frequency 
fi termed the driving frequency in the following. Intro- 
ducing left- as well as right-moving fermions indicated by 
a label L/R and linearizing the dispersion relation around 
the respective Fermi points one arrives at the following 
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This Hamiltonian differs from the equilibrium case only 
through the periodic time dependence of the interac- 
tion U(x — x';t) with an associated Fourier transform 
U q (t) — J dx e~ tqx U(x;t). We assume a repulsive inter- 
action potential U q (t) > of finite range such that it 
is cut off beyond some momentum scale q c . The colons 
: . . . : denote normal ordering relative to the Fermi sea 
and vf is the Fermi velocity. The fermionic density 
Pn(x) with rj — L/R is determined by the fermionic fields 
ipr,(x) = y/2n/LJ2k e~ lkx c kri via p„(x) = : if)l(x)ip n {x) : . 
The operator cjL creates a fermion of the species r? = L/R 
with wave vector k . 

Note that we restrict to the case of spinless fermions. 
This reduced model system already incorporates most of 
the characteristic features of interacting fermions in ID. 
For a system of fermions with spin in equilibrium, for 
example, the dynamics separates into two independent 
sectors of spin and charge, a phenomenon called spin- 
charge separation, each of which can be modeled by a 
Hamiltonian of the form in Eq. (1). Note that the influ- 
ence of a periodic modulation of the Fermi velocity onto 
the momentum distribution in case of fermions with spin 
has been investigated recently 12 . 

In the nonequilibrium scenario under investigation the 
system is initially prepared in the ground state \ipo) at 
some fixed interaction strength such that U q (t < 0) = 
V q (l+v) is chosen to be time- independent for times t < 0. 
At time t — the periodic driving is started with the 
following parametrization of the time-dependence of the 
interaction potential 



U q (t) = V g (1 + v co$(ttt)) . 



(2) 



The dimcnsionless coupling v of the periodic driving is 
chosen v < 1 such that the interaction remains repulsive 
for all times. The interaction strength is characterized 
by the dimensionless number 



2ttvf 
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which within the validity of the Luttinger model is always 
chosen a < 1. We assume that the interaction potential 
V q is cut off beyond the momentum scale q c . In our nu- 
merical simulations we choose a Gaussian for simplicity, 
i.e., V q /{2nv F ) = a exp{- (q / q c ) 2 } . 

Although the Hamiltonian in Eq. (1) is quartic in 
fermionic operators it can be mapped onto a quadratic 
and exactly solvable problem using the bosonization tech- 
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for each right- and left-moving branch r\ = L/R, the 
Hamiltonian in Eq. (1) can be mapped onto a quadratic 
but time-dependent bosonic problem 
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The dispersion of the diagonal part of the above Hamil- 
tonian is given by 



LUq(t) = qV F 

The overall constant A(t) 
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effect on the time evolution of observables except the in- 
ternal energy itself as discussed in Sec. IV. In principle, it 
is possible to diagonalize this Hamiltonian using a time- 
dependent unitary transformation. However, it turns out 
to be suitable to determine the dynamics in the untrans- 
formed basis, see Sec. Ill below. 



III. NONEQUILIBRIUM TIME EVOLUTION 

For the dynamics of all quantities considered such as 
energy density, fermionic density, and the momentum dis- 
tribution it is sufficient to solve the Heisenberg equations 
of motion for the bosonic operators 



dt qv 



(t) 



-iu q (t)b m {t) + iq^r-bl % 



(8) 



with fj the conjugate species of rj, i.e., L = R and vice 
versa. These differential equations for operators can be 
transformed into differential equations for complex func- 
tions Xqrt(t) an d \ri(t) defined by 



bqr,(t) = Xqn{t)bq V + \q V {t)b\ 



(9) 



when inserted into Eq. (8). The resulting system of cou- 
pled differential equations can be cast into a more famil- 
iar form by regarding appropriate superpositions 



>q>) 



Xqv ^qni fiqn — Xqr] H~ A 



qrf 
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The function a qv is the solution of a parametrically 
driven harmonic oscillator and obeys a Mathieu equation 
in properly scaled parameters 



^^+4 [l + 2 7g cos(2r)]a ?J) (r) 
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with the dimensionless time r = Q.t/2, the natural fre- 
quency e q — 2vf1^ 1 \/i + Vq/iftVF) of the harmonic os- 
cillator, and 7 g = isV q /(2Trvp+2V q ) the coupling strength 
of the periodic perturbation. The initial conditions for 
the solution of the Mathieu equation are a qv (t = 0) = 1 
and Oi' qr) {t = 0) = —ivpq- The remaining function 
Pqv = da qr) (t) / dt is proportional to the time 

derivative of a qT] . 

Concluding, the time evolution in the periodically 
driven Luttinger liquid is equivalent to a set of paramctri- 
cally driven harmonic oscillators. The Mathieu equation 
in Eq. (11) in general exhibits no analytic solution in 
terms of elementary functions. For special cases, how- 
ever, such as parametric resonance, approximate analyt- 
ical solutions are available, see below. For the general 
case we solve the differential equations numerically using 
a standard 4-th order Runge-Kutta algorithm. 

The driven harmonic oscillator in Eq. (11) shows an 
instability with exponentially growing amplitudes in the 
case of parametric resonance which occurs for that partic- 
ular plasmonic mode q* for which e q * = 1 or equivalently 

Sl/2 = v F q* ^/l + V r /{nv F ) 18 . 

In the following it is important to distinguish two dif- 
ferent cases of fast and slow driving. The energy scale 
ST associated with the crossover between the two limits 
is set by 



sr 



VFQc 
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For slow driving SI S7* the resonant bosonic mode q* 
for which parametric resonance occurs is determined by 



q* f2<fT 1 



(13) 



to leading order in fi/Q*. Note, however, that not only 
q* but also momenta q within a finite interval of nonzero 
length contribute to the resonance 18 . The rate T of the 
associated exponential growth in time can be determined 
using standard methods 18 



4 1 + 2a 
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The time scale t* for the onset of the instability is then 
determined by the rate T via 



t* 



(15) 



In the opposite case SI 3> SI* of fast driving the resonant 
mode 



q* n»n* SI 
20* 



(16) 



is independent of the interaction potential up to correc- 
tions suppressed by the cutoff q c . The associated rate of 
the exponential growth is then given by 



r 0> ^* lfw„. 



(17) 



Its precise behavior for fi/O* ^> 1 or equivalently 
q* I q c 3> 1 depends on the details of the large momentum 
behavior of the interaction potential. If V q ~ exp[— Cq] 
for some constant C > the rate T ~ exp[— DQ/Cl*] with 
D = Cq c /2 is suppressed exponentially Analogously, 
algebraically decaying potentials V q ~ (<z/9c) _AI yield a 
power law dependence T ~ (S1/ST) 1_A1 for // > 1. 

As r is linear in SI for slow driving, see Eq. (14), the 
adiabatic approximation breaks down after a fixed num- 
ber of periods N pel — fl/T = 4(1 + 2a) /(av) set by the 
interaction strength a and the coupling to the periodic 
perturbation v irrespective of SI. By reducing the driving 
frequency one cannot increase the number of periods for 
the validity of the adiabatic approximation. 



IV. INTERNAL ENERGY DENSITY 

Typically, the periodically driven quantum many-body 
systems considered so far in the literature ignore the pos- 
sible influence of dissipation mechanisms onto the dy- 
namics, with the exceptions of Refs. 8 ' 10 . This is an im- 
portant issue because the energy in the system will in 
general increase during the considered nonequilibrium 
protocol. For bounded Hamiltonians, fermionic or spin 
systems on a lattice, for example, the internal energy 
will necessarily saturate, for unbounded Hamiltonians 
this need not be the case. Even though the internal en- 
ergy for bounded Hamiltonians will always stay finite, 
the question whether the unavoidable presence of dissi- 
pation mechanisms may at some point in time have a 
considerable influence onto the dynamics is still largely 
unanswered. Due to the parametric instability the inter- 
nal energy of the Luttinger liquid diverges as we will show 
below. The existence of the instability naturally sets the 
time scale t* , see Eq. (15), beyond which additional inter- 
nal properties such as the curvature of the fermionic dis- 
persion relation or external dissipation mechanisms have 
to be included for a realistic description. 

The internal energy density £(t) = 
L~~ 1 (ipo(t)\H(t)\ip (t)) of the periodically driven 
Luttinger liquid system at time t is given by 



£{t) 



where 
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9 g ) + |A gL (i)| 2 cosh(20 9 ) 
+ Re[ XqL (t)X* qL (t)] sinh(2^), 
sinh(2^)Re[A2 L+x 2 L ] + 

+2cosh(26» g )Re[A gL XgL], 



(19) 



and (jJ q (t) given by Eq. (7). The Bogoliubov angles 
6 q for the initial state arc determined by the formula 
tanh(26» (? ) = V q (l + v)/{2TTv F + V q {l + v)). Note that the 
Luttinger liquid Hamiltonian in Eq. (1) appears typically 
as the low-energy theory derived from more complicated 
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FIG. 1: (color online) Time evolution of the internal energy 
density £ (t) (blue lines) and its time average £(t) (red dashed 
lines) for different driving frequencies Q in units £* — q c Q* . 
For these numerical simulations we have chosen an analytic 
interaction potential V q / (2-KVp) = aexp(—(q/q c ) 2 ) with a = 
1/2 and v — 1/5. The zero of energy has been chosen such 
that £ (t = 0) = 0. 



energy photons 6 . In the present model system this mech- 
anism is provided due to the finite range q c of the interac- 
tion. Physically speaking, the system is not able to follow 
the fast external perturbation and only perceives its aver- 
age contribution. As the initial state is not an eigenstate 
of the time-averaged Hamiltonian for times t > the 
dynamics becomes equivalent to that of an interaction 
quench. Indeed, we find that the time-averaged energy 
density £ (t) follows precisely the behavior of the interac- 
tion quench scenario. 

When the driving frequency is lowered the rate of the 
instability T grows and thus the onset of the instability 
moves to smaller times t* . For fl ss fi* the transient dy- 
namics is directly followed by an exponentially increasing 
contribution due to the influence of the instability. In this 
way the instability is so strong that it hinders the buildup 
of a meta-stable state completely. For even lower frequen- 
cies, however, the rate T ~ f2 decreases again opening up 
the window for the meta-stable state as can be clearly 
observed in Fig. 1. 



V. FERMIONIC DENSITY 



many-body systems. Within such a mapping additional 
time-dependent contributions to A(t) in Eq. (6) can be 
generated in the periodically driven case. Those contri- 
butions depend on the details of the model and have to 
be worked out for each particular case. In the present 
work we are interested in the generic low-energy prop- 
erties that are all contained in the Luttinger liquid de- 
scription of the model Hamiltonian in Eq. (1). Thereby 
we ignore additional contributions to the energy density 
that are generated by the mapping onto this low-energy 
theory. 

In Fig. 1 plots for £(t) are shown for different driving 
frequencies. Periodically driven systems without insta- 
bilities develop stationary states at long times where ex- 
pectation values time-averaged over one period become 
time-independent. Besides the internal energy density we 
have included its time average 

I(t) = —j dt'£{t') (20) 

in Fig. 1 for the identification of such stationary states. 

In case of fast driving fi 3> the rate T of the para- 
metric resonance is strongly suppressed due to the finite 
cutoff scale q c of the interaction potential, see Eq. (17). 
After the initial transient dynamics following the start 
of the periodic driving the system settles to a stationary 
state with a constant time-averaged energy density. As 
r is small but still finite the system will nevertheless de- 
velop the instability for times t > t* . Consequently this 
intermediate state is only meta-stable. 

In the limit f2 — > oo the dynamics becomes effectively 
equivalent to one of a time-averaged Hamiltonian if there 
exists a mechanism that prevents the absorption of high- 



The internal energy density of the system mirrors the 
instability by showing a divergence at times t > t* . The 
question, however, which internal perturbations can pre- 
vent the buildup of the instability we have not touched 
up to now. Regarding the time evolution of the fermionic 
density we argue that it may be sufficient to include the 
curvature of the fermionic dispersion relation. 

In the following we analyze the dynamics of an initially 
localized fermionic density wave packet in presence of the 
periodic driving. The time evolution of a local perturba- 
tion in the fermionic density is solely determined by the 
solution a qrj (t) of the Mathieu equation, see Eq. (11), 

p(x, t) = 2j ^ cos{qx) p° q Re[a qL (t)] (21) 

with Pq = J dqe~ lqx po{x) characterising the initial den- 
sity profile po(x) — p(x,t = 0). In Fig. 2 we show for 
the fast and slow driving regimes the dynamics of the 
fermionic density where we have chosen a Gaussian wave 
packet as initial condition for illustration. 

For fast driving Q 3> fi* the initial local perturbation 
splits into right- and left-moving contributions. As in 
case of the internal energy density in Sec. V the dynamics 
becomes equivalent to the interaction quench scenario. 
This picture is suitable for times t < t* before the onset 
of the instability. In this regime the time scale t* oc f2 _1 , 
see Eq. (14), is large and grows in a power-law fashion 
for small values of the driving frequency. 

In the opposite case fl <C f2* a completely differ- 
ent picture emerges. After a fixed number of periods 
AT per ~ n/r = 4(1 + 2a)/{av) to leading order in fl/n* 
the dynamics is dominated by the parametric resonance 
leading to an exponential growth of the initial perturba- 
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FIG. 2: (Color online) Time evolution a localized perturba- 
tion in the fermionic density p(x, t) for the cases of fast (upper 
plot) and slow (lower plot) driving. Again, we have chosen an 
analytic interaction potential V q /(2nvF) = aexp(—(q/q c ) 2 ) 
with a = 1/2 and v = 1/5 for these numerical simula- 
tions. The initial density profile is a Gaussian p(x, t = 0) = 
po exp(—(x/x c ) 2 )/2 with x c q c = 2 in the slowly and x c q c = 10 
in the fast driven case. The buildup of the superlattice in the 
slowly driven limit is independent of the precise choice of the 
width x c of the initial density wave packet. In the fast driv- 
ing case, that is equivalent to an interaction quench, see main 
text, the universal low-energy limit corresponds to x c q c 3> 1. 

tion for x <C v max t 

p(x, t) q c A cos{q*x) cos + e Tt (22) 

with A a constant nonuniversal prefactor. A superlattice 
forms with a period 

x* = 2n/q* (23) 

set by the resonant mode, see Eqs. (13,16), whose am- 
plitude is growing exponentially at a rate T. Note that 
this does not lead to a violation of particle number con- 
servation as it might seem from the plot in Fig. 3. For 
each density hump there also exists a valley of depletion 
of fermionic charge carriers. Moreover, the superlattice 
extends only over a distance d ~ v max t within the light 



cone set by the maximal sound velocity. Integrating over 
the whole real-space shows that the particle number is 
still conserved as one can directly check via Eq. (21). 

Now we want to argue that including the curvature 
of the fermionic dispersion relation will cut off the ex- 
ponential growth of the superlattice for sufficiently large 
densities. In the limit of slow driving we can approx- 
imately neglect the influence of the finite range of the 
interaction and set V q ~ Vq. For a q-independent in- 
teraction the influence of a nonlinear fermionic disper- 
sion relation can be accounted for approximately 19 . The 
time-independent version of the Hamiltonian in Eq. (1) 
including the quadratic curvature contribution can be 
mapped to a free Fermi gas 19 . For a free Fermi gas the 
buildup of large densities is prevented due to the non- 
linear dispersion eventually leading to the production of 
shock waves 20 . 



VI. MOMENTUM DISTRIBUTION 

In equilibrium the momentum distribution for the 
fermionic particles exemplifies the different influence of 
repulsive interactions in one dimension compared to 
higher dimensions where Fermi liquid theory holds. In 
Luttinger liquids the momentum distribution shows no 
jump at the Fermi energy even at zero temperature rem- 
iniscent of the absence of a finite quasiparticle weight. 

The momentum distribution function fk (t) for the left- 
moving fermions 

/*(*) - (Mc\ L {t)c kL (t)\^), (24) 

with l^o) the initial state, is connected to an equal- 
time correlation function f{x,t) in real-space via Fourier 
transformation 

/*(*) = / |V te /0M), (25) 

f(x,t) = (i> \Tp[(x,t)i> L (0,t)\i> ) (26) 

that can be calculated analytically using the bosonization 
technique 17 

/(*,*) = — ^- e-*W>. (27) 
a + ix 

with a -1 an ultraviolett cut-off. The rate function 
F(x,t) = Feq(x) + F p (x,t) can be separated into an equi- 
librium part F eq associated with the initial state and a 
nonequilibrum contribution F p (x,t) due to the periodic 
driving. 

F„(x) = 4jT^sm 2 (f) S mh 2 (fy, (28) 

FM-<£ |^(f) 

x [\X qL \ 2 cosh(20„) + Re( XgL A; L ) sinh(20 g )] 
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FIG. 3: (Color online) Time evolution of the momentum dis- 
tribution fk(t) for fast (upper plot) and slow (lower plot) driv- 
ing. The parameters of the numerical simulation are chosen 
as for the previous plots. 



Here, 9 q denotes the Bogoliubov angles of the diagonal- 
izing transformation for the initial equilibrium Hamilto- 
nian obeying the equation tanh(2# g ) = V q (l+v) / (2ttvf + 
V q (l + v)). In Fig. 3 we show numerical results for the 
momentum distribution for fast and slow driving. 

For Q 3> £1* the behavior under time evolution is con- 
sistent with the picture observed for the internal energy 
density in Sec. IV. Following the initial transient dynam- 
ics the momentum distribution becomes meta-stable on 
intermediate times t < t* . The momentum distribution 
is self-averaging in the sense that it is time-independent 
in contrast to the energy density where only the time 
average becomes constant. In the limit £1 — > oo we re- 
cover the interaction quench limit of an effectively time- 
averaged Hamiltonian as for the internal energy density, 
see Sec. IV. In the vicinity of the Fermi level the momen- 
tum distribution shows a nonanalytic behavior for large 
times 



rik(t —> oo) 



1 n- 
2 




FIG. 4: (Color online) Time evolution of the fermionic de- 
fect density nd c f(i) (blue line) and its time average rTd c f(t) 
(red dashed line) over one period. The parameters for the 
numerical simulations have been chosen as before. 



with sgn(fc) the sign function. The exponent 7 is in pre- 
cise agreement with an interaction quench scenario 21 . 

For slow driving O 17* the momentum distribution 
develops steps as can be seen in Fig. 3, see also Ref. 13 . A 
similar observation has been made recently for the case 
of a periodic modulation of the Fermi velocity in a Lut- 
tinger liquid 12 . These steps may be associated with scat- 
tering processes between fermions under the absorption 
of quanta of the driving frequency fl. Thus, the dominant 
processes under the periodic driving are not only energy- 
conserving ones but also those where energy is conserved 
up to multiples of the driving frequency. Note that this 
step structure is remarkably similar to a simplified pic- 
ture where in spirit of the work by Tien and Gordon 22 for 
noninteracting systems the periodic driving generates a 
weighted superposition of equilibrium momentum distri- 
butions shifted by an energy nfi. It is, however, not pos- 
sible to establish such a superposition principle precisely 
in the present interacting system. Due to the paramet- 
ric instability the fermions are redistributed completely 
for times t — ¥ 00 leading to a momentum distribution 
fk(t —> 00) = 1/2 of infinite temperature. 

Due to scattering processes induced by the periodic 
driving fermions are redistributed in the vicinity of the 
Fermi surface already for small times. This excludes the 
possibility for an adiabatic limit where the momentum 
distribution follows adiabatically the external perturba- 
tion. This generates fermionic defects 



n de f(t) = 



00 dk 
2tt 



(30) 



sga(k)\k/q c \ 



(29) 



compared to the adiabatic limit. Here, f^(t) is the equi- 
librium momentum distribution for the instantaneous in- 
teraction potential U q (t). The integral can be restricted 
to positive momenta as fk(t) — fjl q (t) = — [/-fc(t) — 
/-*!(*)]• A s shown in Fig. 4 after some transient dy- 
namics the time-averaged defect density n dc f(t) increases 
linearly in a regime where in contrast the time-averaged 
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energy density £ (f) settles to a constant value on inter- 
mediate time scales which we have identified as a meta- 
stable state. Even though £ (t) approaches a constant 
value interaction energy is successively transferred to ki- 
netic energy reminiscent of a heating process with an 
increasing entropy. 

VII. CONCLUSIONS 

In this work we have studied the dynamics of a Lut- 
tinger liquid with a periodically time-dependent repul- 
sive interaction potential. Under the periodic driving the 
system develops an instability due to a parametric res- 
onance. The associated time scale t*, see Eqs. (14,17), 
sets the limit beyond which dissipation mechansims or 
internal perturbations have to be included into an ap- 
propriate description of the dynamics. On intermediate 
time scales before the onset of instability, it is possible to 



identify meta-stable states for fast and slow driving with 
constant time-averaged internal energies. The paramet- 
ric instability generates an exponential growth of pertur- 
bations in the fcrmionic density leading to the buildup 
of a superlattice with period x* , confirm Eq. (23). The 
fermionic momentum distribution develops a step struc- 
ture that can be associated with fermiomic scattering 
processes under the absorption or emission of quanta of 
the driving frequency f2. 

Acknowledgments 

We thank S. Kehrein for fruitful discussions. This work 
was supported by the Center for Nanoscience CeNS Mu- 
nich and the Deutsche Forschungsgemeinschaft through 
SFB 631 and SFB TR 12. Financial support by the Ex- 
cellence Cluster "Nanosystems Initiative Munich (NIM)" 
is gratefully acknowledged. 



1 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, 
D. Abusch-Magder, U. Meirav and M. A. Kastner, Nature 
391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp and 
L. P. Kouwenhoven, Science 281, 540 (1998); J. Schmid, 
J. Weis, K. Eberl and K. von Klitzing, Physica B 258, 
182 (1998): W. G. van der Wiel, S. De Franceschi, T. Fuji- 
sawa, J. M. Elzerman, S. Tarucha and L. P. Kouwenhoven, 
Science 289, 2105 (2000). 

2 I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

3 A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalat- 
tore, Rev. Mod. Phys. 83, 863 (2011). 

4 P. Nordlander, N. S. Wingreen, Y. Meir, and D. C. Lan- 
greth, Phys. Rev. B 61, 2146 (2000); R. Lopez, R. Aguado, 
G. Platero, and C. Tejedor, Phys. Rev. Lett. 81, 4688 
(1998); T. K. Ng, Phys. Rev. Lett. 76, 487 (1996); 
M. H. Hettler and H. Schoeller, Phys. Rev. Lett. 74, 4907 
(1995). 

5 A. Kaminski, Y. V. Nazarov, and L. I. Glazman, Phys. 
Rev. Lett. 83, 384 (1999); A. Kaminski, Y. V. Nazarov, 
and L. I. Glazman, Phys. Rev. B 62, 8154 (2000); 
Y. Goldin and Y. Avishai, Phys. Rev. Lett. 81, 5394 
(1998); Y. Goldin and Y. Avishai, Phys. Rev. B 61, 16750 
(2000)1; A. Schiller and S. Hershfield, Phys. Rev. Lett. 77, 
1821 (1996). 

8 M. Heyl and S. Kehrein, Phys. Rev. B 81, 144301 (2010). 

7 J. K. Freericks, V. M. Turkowski, and V. Zlatic, Phys. Rev. 
Lett. 97, 266408 (2006); J. K. Freericks, Phys. Rev. B 77, 
075109 (2008); N. Tsuji, T. Oka, and H. Aoki, Phys. Rev. 
B 78, 235124 (2008). 

8 N. Tsuji, T. Oka, and H. Aoki, Phys. Rev. Lett. 103, 
047403 (2009); 

9 A. V. Joura, J. K. Freericks, and Th. Pruschke, Phys. 
Rev. Lett. 101, 196401 (2008); M. Eckstein and P. Werner, 



Phys. Rev. B 82, 115115(2010); M. Eckstein, T. Oka, and 
P. Werner, Phys. Rev. Lett. 105, 146404(2010); N. Tsuji, 
T. Oka, P. Werner, and H. Aoki, Phys. Rev. Lett. 106, 
236401 (2011); M. Eckstein and P. Werner, Phys. Rev. 
Lett. 107, 186406 (2011); 

10 C. Aron, G. Kotliar, and C. Weber, Phys. Rev. Lett. 108, 
086401 (2012); A. Amaricci, C. Weber, M. Capone, and 
G. Kotliar, arXiv:1106.3483v3 (2012). 

11 T. Oka and H. Aoki, Phys. Rev. Lett. 95, 137601 (2005); 
T. Oka and H. Aoki, Phys. Rev. B 78, 241104(R) 

(2008) ; M. Mierzejewski and P. Prelovsek, Phys. Rev. 
Lett. 105, 186405 (2010); M. Mierzejewski, J. Bonca, and 
P. Prelovsek, Phys. Rev. Lett. 107, 126601 (2011). 

12 C. D. Graf, G. Weick, and E. Mariani, Europhys. Lett. 89, 
40005 (2010). 

13 Yu. Kagan and L. A. Manakova, Phys. Rev. A 80, 023625 

(2009) . 

14 S. Pielawa, Phys. Rev. A 83, 013628 (2011). 

15 A. Eckardt, C. Weiss, and M. Holthaus, Phys. Rev. Lett. 
95, 260404 (2005). 

16 D. Poletti and C. Kollath, Phys. Rev. A 84, 013615 (2011). 

17 J. von Delft and H. Schoeller, Ann. Phys. (Leipzig) 7, 225 
(1998). 

18 L. D. Landau and E. M. Lifshitz, Mechanics (Elsevier, Am- 
sterdam, 2008). 

19 A. V. Rozhkov, Eur. Phys. J. B 47, 193 (2005); A. Imam- 
bekov and L. I. Glazman, Science 323, 228 (2009). 

20 E. Bettelheim, A. G. Abanov, and P. B. Wiegmann, Phys. 
Rev. Lett. 97, 246402 (2006). 

21 M. A. Cazalilla, Phys. Rev. Lett. 97, 156403 (2006); 
J. Rentrop, D. Schuricht and V. Meden, New J. Phys. 14, 
075001 (2012). 

22 P. K. Tien and J. P. Gordon, Phys. Rev. 129, 647 (1963). 



